Real-Time Parameter Estimation in the 
Frequency Domain 

E. A. Morelli 


Reprinted from 


Volume 23, Number 5, Pages 812-818 



A publication of the 

American Institute of Aeronautics and Astronautics, Inc. 
1801 Alexander Bell Drive, Suite 500 
Reston, VA 20191-4344 



Journal of Guidance, Control, and Dynamic s 
Vol. 23, No. 5, September-October 2000 


Real-Time Parameter Estimation in the Frequency Domain 

Eugene A. Morelii* 

NASA Langley Research Center, Hampton, Virginia 23681-2199 


A method for real-time estimation of parameters in a linear dynamic state-space mode! was developed and 
studied. The application is aircraft dynamic model parameter estimation from measured data in flight Equation 
error in the frequency domain was used with a recursive Fourier transform for the real-time data analysis. Linear 
and nonlinear simulation examples and flight test data from the F-18 High Alpha Research Vehicle were used to 
demonstrate that the technique produces accurate model parameter estimates with appropriate error bounds. Pa- 
rameter estimates converged in less than 1 cycle of the dominant dynamic mode, using no a priori information, with 
control surface inputs measured in flight during ordinary piloted maneuvers. The real-time parameter estimation 
method has low computational requirements and could be implemented aboard an aircraft in real time. 


Nomenclature 

A, B,C,D = system matrices 

E{ } = expectation operator 

g = acceleration due to gravity, ft/s 2 

h = altitude, ft 

j = imaginary number, ,/(—!) 

M = Mach number 

N = total number of samples 

R - measurement noise covariance matrix 

Re = real part 

t = time 

V t = true airspeed, ft/s 

x , «,y = state, input, and output vectors 

Zi = measured output vector al time /' A i 

a - angle of attack, rad 

p = sideslip angle, rad 

S u , & r = aileron, rudder deflections, rad 

<5,, 8, = elevator, stabilator deflections, rad 

Sjj = Kronecker delta 

9 = p-dimensional parameter vector 

u, = discrete measurement noise vector 

or 2 = variance 

10 = angular frequency, rad/s 

Subscripts 

i ss value at time / At 

0 = trim or initial value 

Superscripts 

T - transpose 

t = complex conjugate transpose 

= discrete Fourier transform 
* = estimate 

— I = matrix inverse 

Introduction 

R EAL-TIME identification of dynamic models is a requirement 
for indirect adaptive or reconfigurable control. 1 One approach 
for satisfying this requirement is to assume the dynamic model has a 
linear structure with time- vary ing parameters to account for changes 


Received 10 August 1999; revision received 1 January 2000; accepted for 
publication 13 February 2000. Copyright © 2000 by the American Institute 
of Aeronautics and Astronautics, Inc. No copyright is asserted in the United 
States under Title 17, U S Code The U S. Government has a royalty-free 
license to exercise all rights under the copyright claimed herein for Govern- 
mental purposes. All other rights are reserved by the copyright owner. 

* Research Engineer, Dynamics and Control Branch. 


in the flight condition, slores, configuration, remaining fuel, or from 
various types of failures, wear, or damage. The task is then to identify 
accurate linear model parameter estimates from measured data in 
real time, so that the adaptive control logic can make the necessary 
changes to the control law to achieve stability and performance 
goals. 

Two main problems plague accurate real-time parameter estima- 
tion: noise and data information content. It is difficult to design a 
parameter estimation technique that is insensitive to noise but still 
responds rapidly to sudden changes in the system dynamics, mainly 
because it takes a fairly long data record to distinguish noise from a 
sudden change in the dynamics. This problem can be handled in the 
time domain using recursive least squares and a forgetting factor, 2 
or by using sequential batch least squares with short data records 
and including various constraints in the parameter estimation cost 
function. 3 4 If an extended Kalman filtering approach 5 8 is used, 
discriminating signal from noise is implemented through weighting 
matrices that represent assumed measurement and process noise co- 
variances. For all of these time-domain methods, some adjustment 
of one or more tuning parameters must be done in simulation. In ad- 
dition, the standard errors for the model parameter estimates, which 
are important both for failure detection and adaptive or reconfig- 
urable control, cannot be accurately and reliably computed using 
recursive or sequential batch time-domain methods. 

In the context of airplane flight, lack of information content in 
the data can be problematic because there are frequently extended 
periods where the control and state variables are fairly constant. 
Signal levels are at or below the (relatively constant) noise level. In 
this circumstance, a time-domain regression method will give very 
inaccurate parameter estimates unless the estimation is regularized 
by including a term in the cost function that penalizes movement of 
the parameters away from some a priori known values (for example, 
values from wi nd-tunnel tests) and/or a term that penali zes ti me vari- 
ation of the parameter estimates. 4 Tuning parameters are required 
for this approach because the magnitude of the penalty term(s) must 
be balanced properly relative to the least-squares part of the cost 
function used for parameter estimation based on measured data. 

Another problem that falls in the category of poor data informa- 
tion content is data collinearity due to the control system. 9 Many 
control laws move more than one control surface at the same time or 
move control surfaces in proportion to state variables with a small 
time delay. When states and controls are nearly proportional to one 
another, it is impossible to identify individual stability and control 
derivatives from the measured data alone. 

There are many parameter estimation methods, but the require- 
ment of being simple enough to be implemented in real time aboard 
the aircraft narrows the field. In particular, any method that iterates 
through the data must be eliminated. The current work is an inves- 
tigation of a single-step frequency-domain method for the real-time 
parameter estimation task and an evaluation of its suitability for 
aircraft problems. This real-time parameter estimation method was 
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first proposed as a component of a technique for in-flight system 
identification. 10 In the present work, the real-time parameter estima- 
tion method is developed further and applied to realistic simulation 
and flight-test data. 

The next section gives the problem statement and outlines the 
necessary theory. Following this, the real-time parameter estima- 
tion method is applied to a simulation example, where a linear truth 
model is used with outputs corrupted by noise similar to that seen in 
flight and by noise that is worse than usual. The application is iden- 
tifying an accurate model for the longitudinal rigid-body dynamics 
of a conventional fighter. The real-time parameter estimation pro- 
cedure is further demonstrated with a longitudinal example using 
flight-test data and a nonlinear simulation example, where both lon- 
gitudinal and lateral/directional parameters are estimated. 

Theoretical Development 

Airplane dynamics can be described by the following linear model 
equations 11 : 


x(t) — Ax{t) + Bu(t) (1) 

x(0)~x 0 (2) 

y(t) =Cx(t)+Du(t) (3) 

1 = 1, 2,..., A (4) 


Matrices A, B, C, and D in Eqs. (1) and (3) contain stability and 
control derivatives that are assumed as constant model parameters 
to be estimated from flight data. Repeating the parameter estimation 
at short intervals produces piecewise constant estimates for time- 
varying model parameters in the linear model structure. The input 
quantities are control surface deflections, <5, or 8 S , <5 a , and <5,, with 
the states selected from air data, V',, or, and body axis angular 
velocities, p , q, and r, and Euler angles, 0, 6 , and 0. Output quan- 
tities can include the states and translational accelerations, a x , a y , 
and a z . 

Equation Error in the Frequency Domain 

The finite Founer transform of a signal x{t) is defined by 


an equation error formulation to estimate the stability and control 
derivatives contained in matrices A, B, C, and D . 

For the /;th state equation of vector Eq. (9), the cost function is 

h = ^ Y'. I. i‘»nh(n) - A k x(n) - B k u(n)\ 2 (1 1) 

n = 1 


where A* and B k are the kih rows of matrices A and B , respectively, 
and x k («) is the klh element of vector Jr for frequency co n . Symbols 
x (n) and u(n) denote the Fourier transform of the state and control 
vectors for frequency to n . There are m terms in the summation, corre- 
sponding to m frequencies of interest, and each transformed variable 
depends on frequency. Similar cost expressions can be written for 
individual output equations from vector Eq. (10). 

Denoting the vector of unknown model parameters in A* and B k 
by 6 . the problem can be formulated as a standard least-squares 
regression problem with complex data, 

Y = X0 + e (12) 


where 


joj ix k (\) 
j(D 2 ** (2) 

jv m Xk (w) 


(13) 


i r (l) « r (D 

x r (2) u T ( 2) 

. (14) 

Jc 7 (in) u T (m)_ 

and e represents the complex equation error in the frequency domain. 
The least-squares cost function is 

J = ~ X6)\Y - X0) (15) 


,?(<w) 


f 


x(t)e~ ju)l d/ 


which can be approximated by 


which is identical to the cost in Eq. (11). The parameter vector 
(^) estimate that minimizes this cost function is computed from 13 

0 = [RelA^A') ]~ 1 Rc(A" + F) (16) 


n - i 

i(£«) S3 Ar V Xj e J “'‘ 
/ — 0 


( 6 ) 


The estimated parameter covariance matrix is 

cov(0) = E{(0 - 6)(0 - G ) T ) = or 2 [Re(^r t J0]‘ l (17) 


Subscript i indicates the variable value at lime iAi y and At is the 
sampling interval. The summation in Eq. (6) is defined as the discrete 
Fourier transform, 

v - 1 

X(o» = (7) 

i=0 

so that 


where the equation error variance a 2 can be estimated from the 
residuals, 

a 2 = U/(m - p)][(Y - XG)\Y - X0)\ (18) 

and p is the number of elements in parameter vector 9 Parameter 
standard errors are computed as the square root of the diagonal 
elements of the cov(0) matrix from Eq. (17) using & 2 from Eq (18). 


x(co) ~ A (o>)Af (8) 

Some fairly straight forwaid corrections 12 can be made to Eq. (8) 
to remove the inaccuracy resulting from Eq. (8) being a simple Euler 
approximation to the finite Founer transform of Eq. (5). However, 
if the sampling rate is much higher than the frequencies of interest 
(as is true in this case), then the corrections are small and can be 
safely ignored. 

Applying the Founer transform to Eqs. (1) and (3) gives 


jcox(co) — Ax(co) H- Bu {to) 

(9) 

y(a>) — Cx(ari + Du{a>) 

(10) 


When the states, outputs, and inputs are measured, individual 
state or output equations from vector Eqs. (9) or (10) can be used in 


Recursive Fourier Transform 

For a given frequency to, the discrete Fourier transform in Eq. (7) 
at sample time / is related to the discrete Fourier transform at time 
/ — 1 by 

Xi(w) — Xj _ j {to) -F Xj exp(— jwiAt) (19) 

w'hcre 

cxp(- jcoiAt) = cxp(— jcoAt) exp — 1) A/] (20) 

The quantity c\p{--jcoAt) is constant for a given frequency and 
constant sampling interval, it fol low’s that the discrete Fourier trans- 
form can be computed for a given frequency at each time step using 
one addition in Eq. (19) and two multiplications: one in Eq. (20) 
using the stored constant exp(- juoAt ) for frequency to, and one 
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in Eq. (19) There is no need to store the time-domain data in 
memory when computing the discrete Fourier transform in this way 
because each sampled data point is processed immediately. Time- 
domain data from all preceding maneuvers can be used in all subse- 
quent analysis by simply continuing the recursive calculation of the 
Fourier transform. In this sense, the recursive Fourier transform acts 
as memory for the information in the data. More data from more 
maneuvers improves the quality of the data in the frequency domain 
without increasing memory requirements to store it. In addition, the 
Fourier transform is available at any time /A/, The approximation 
to the finite Fourier transform is completed using Eq. (8). 

Rigid-body dynamics of piloted aircraft lie in the rather narrow 
frequency band of approximately 0.01—1.5 Hz. Theretore, it is pos- 
sible to select closely spaced fixed frequencies for the Fourier trans- 
form and the subsequent data analysis. Excluding zero trequency 
removes trim values and measurement biases, and so it is not nec- 
essary to estimate bias parameters. Using a limited frequency band 
for the Fourier transformation confines the data analysis to the fre- 
quency band where the system dynamics reside and automatically 
tillers wideband measurement noise or structural response outside 
the frequency band of interest. 

f or ah of the examples studied in this work, frequency spac- 
ing of 0.04 Hz on the interval [0 1-1.5] Hz was found to be ad- 
equate, giving 36 evenly spaced frequencies lor each transformed 
time-domain signal Finer frequency spacing requires slightly more 
computations, but was found to have little effect on the example re- 
sults. When the frequency spacing is very coarse, there is a danger 
of omitting important frequency components, and this can lead to 
inaccurate parameter estimates. In general, a good rule of thumb 
is to use frequencies evenly spaced at 0.04 Hz over the bandwidth 
for the dynamic system. For good results, the bandwidth should be 
limited to the frequency range where the signal components in the 
frequency domain are at least twice the amplitude of the wideband 
noise level m the frequency domain However, the algorithm is ro- 
bust relative to these design choices, so that the selections to be 
made aic not difficult. 

For airpkme dynamic modeling, the number ot time-domain sig- 
nals So be transformed is usually low (nine or less, more if there 
arc many control surfaces so that this approach requires a small 
amount of computei memory. Because the data analysis is done in 
the frequency domain, the memory required is fixed and independent 
of the time length of the flight maneuvers. 

The lecursive Fourier transform update need not be done for 
every sampled time point. Skipping some time points effectively 
decimates the data prior to Fourier transformation. This saves com- 
putation and does not adversely impact the frequency-domain data 
because the Nvquisl frequency (equal to one-half the sampling fre- 
quency^ is usually much higher than the relatively low frequencies 
being used in the recursive Fourier transform. 

Examples 

For longitudinal aircraft short-period dynamics, the state vector 
x and input vector u in Eq. (1) are defined by 

x = \a q\ T u = [8,] (21) 

System matrices containing the model parameters are 



This modei assumes a effects can be subsumed into the Z' and 
M ({ derivatives. Parameter Z[ includes the inertial term, that is, 
Z;' — 1 -p Z, r In this and all of the other examples, state equations 
were used for the equation error parameter estimation. 

In the first example, a perturbation elevator input was applied to a 
know-n linear model to produce simulated state and output responses. 
Figure 1 shows the elevator input <5,. The first 8 s of the elevator 
input w'ere taken from measured flight test data for a longitudinal 
tracking task. Constant elevator deflection is held for the final 7 s, 
to simulate low data information content. The simulated aircraft 
is a conventional F-16 (Refs. 1 1 and 14) with forward c.g. position 
(0.2c) m straight and ievel trimmed flight at 10,000 ft, 7-deg angle of 


attack, and Mach 0.37. The simulated outputs were corrupted with 
20% Gaussian random white noise. This made the signal-to-noise 
ratio 5-to-l for each simulated output measurement. Figure 2 shows 
the simulated measured perturbation angle of attack a and pitch rate 
q. The elevator input was assumed to be measured without noise, 
which is a close approximation to reality. Model parameter values 
used to generate the simulated test data, called true values, are given 
in column 2 of Table I. Parameter estimation was done in real time 
using equation error in the frequency domain applied to the tw'o 
state equations, with the Fourier transform computed recursively, as 
described earlier. All angular quantities were expressed in radians 
for the data analysis, but were plotted in degrees. Angular rates were 
expressed in radians per second for the data analysis and plotted in 
degrees per second. 

Figure 3 shows a time history of M a model parameter esti- 
mates based on the simulated noisy data only. Plots for the other 
model parameters were similar. The model parameter estimates and 


Table 1 F-16 linear simulation results, qo = 7 deg, 

Ao = 10,000 ft, Mq - 0.37 


20% noise estimate. 

Parameter True value (standard error) 

50% noise estimate, 
(standard error) 

Zu 

-0.600 

-0.624 

-0.602 



(0.047) 

(0.132) 


0.950 

0.960 

0.986 



(0.016) 

(0.044) 

7.6, 

-0.115 

-0.104 

-0.134 



(0.017) 

(0.049) 

M u 

-4.300 

-4.197 

-4.021 



(0.136) 

(0.354) 

M (t 

-1.200 

-1.238 

-1.200 



(0.045) 

(0.119) 

Mi, 

-5.157 

-5.157 

-5.246 



(0.048) 

(0.130) 



Time (s) 


Fig. 1 Elevator input. 




Time (s) 

Fig. 2 Simulated measured outputs, 20% noise. 
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Fig. 3 Parameter estimation, 20% noise. 




0 5 10 15 

Time (s) 

Fig. 4 Simulated measured output, 50% noise. 

standard errors were computed once a second using Eqs. 06-18) 
with frequency-domain data from the recursive Fourier transform tn 
Eqs. ( 1 9) and (20). This update rate for the parameter estimation was 
used throughout the examples, but was chosen aibitrarily. Parameter 
estimation updates can be done at a faster or slower rate, w ith the up- 
per limit defined by the rate used in the recursive Fourier transform. 
The algorithm required no starting values for the parameters, and 
the estimates were not regularized in any way with a prion values 
or constraints on temporal changes in the parameter estimates. The 
first parameter estimates are shown at 2 s, because the parameter 
estimates after 1 s were poor with very large standard errors, due 
to the lack of information content in the data during the first sec- 
ond of the maneuver. Initial parameter estimate accuracy and speed 
of convergence could be improved using a prion information, but 
this was avoided so that the performance of the real-time parameter 
estimation algorithm alone could be studied. 

Figure 3 shows that the parameter estimates converge to the true 
value. The calculated standard errors are representative of the es- 
timated parameter accuracy throughout the maneuver and do not 
suffer from the covariance wind-up problem characteristic of recur- 
sive time-domain methods. Column 3 of Table l contains parameter 
estimates with corresponding standard errors in parentheses. These 
results are for the end of the 15-s maneuver. Every parameter esti- 
mate is within ±1 standard error of the true value, indicating that 
the parameter estimation is accurate and that the estimated standard 
errors properly represent the true accuracy of the parameter esti- 
mates. The calculated standard errors accurately convey information 
on the quality of the parameter estimates throughout the maneuver 
and do not become smaller with increasing maneuver time when 
there is no information in the data. 

Figure 4 shows the same simulated model outputs using the same 
input to the same model, but with the added Gaussian random white 
noise level raised from 20% to 50%, reducing the signal-to-noise 
ratio to 2-to- 1 . In addition, two simulated data dropouts with values 
of - 100 were added to the simulated pitch rate output. Plots in Fig. 5 
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Fig. 5 Parameter estimation, 50% noise. 


show that the parameter estimates in this case again converge to the 
true values, although the standard error values are generally higher, 
due to the increased noise level. Plots for the model parameters not 
shown were similar. The convergence rate of the parameter esti- 
mates to the true values was similar to the lower noise case, requir- 
ing approximately 4 s of data This corresponds to approximately 
1 .4 periods of the short-period natural frequency for the simulation 
model. Considering that no substantial information is contained in 
the data for the first 2 s, this is an excellent result. Discounting the 
first 2 s, the parameters were accurately estimated from approxi- 
mately 2 s of data, corresponding to 0.7 cycles of the short-period 
mode. This may not be fast enough for effective use in an adaptive 
or reeonfigurable control scheme, in which case some mechanism 
for augmenting the data information, such as a priori parameter es- 
timates or auxiliary signal injection on the control surfaces, wall be 
required to reduce the convergence time. Still, considering that algo- 
rithm used only simulated measured data from an ordinary piloted 
maneuver, the results show' a fast and accurate convergence for the 
parameter estimates. 

The real-time parameter estimation algorithm is robust to mea- 
surement noise levels and infrequent data dropouts because of the 
automatic filtering inherent in using a limited bandwidth for the 
recursive Fourier transform. In effect, the data dropouts look like 
high-frequency noise. Column 4 of Table 1 gives the parameter es- 
timates and standard errors for the 50%' noise case at the end of 
the 15-s maneuver. As before, every parameter estimate is within 
±1 standard error of the true value, indicating that the parameter 
estimates and standard errors are accurate. The algorithm gave good 
results for even higher noise levels, with the upper limit noise level 
corresponding to a signal-to-noise ratio around 0.5. As noise level in- 
creased further, parameter estimate accuracy gradually deteriorated. 

The linear simulation and the real-time data analysis were pro- 
grammed and run in MATLAB® 5.3 (Ref. 15). Sampling rate for the 
data was 40 Hz, and the recursive Fourier transform updates were 
done at 40 Hz The real-time estimation algorithm ran roughly 10 
times faster than real time (1.5 s for a 15-s maneuver) on a Gate- 
way 450-MHz E-4200 serial processor, running Microsoft Windows 
NT 4.0. 
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The next example used real flight-test data from the F- 1 8 High 
Alpha Research Vehicle (HARV) to demonstrate the real-time pa- 
rameter estimation method. Figure 6 shows the measured stabilator 
deflection for this 14-s maneuver. Measured outputs are shown in 
Fig. 7. The maneuver was a sequence of doublets implemented by 
the pilot from a steady trim condition at 20-deg angle of attack, 
24, 1 00- ft altitude, and Mach 0.34. Parameter estimates are plotted 
as a function of lime in Fig. 8, along with dashed lines indicat- 


F-18 HARV Maneuver 153B 



Fig. 6 Stabilator input. 


F-18 HARV Maneuver 153B 



Fig. 7 Measured outputs. 


ing the ±2 a error bounds. Parameter estimates from standard batch 
time-domain output error parameter estimation 16 are plotted as solid 
lines for comparison. The real-time parameter estimation algorithm 
produced parameter estimates and standard errors that were in agree- 
ment with the batch time-domain estimates after about 5 s, of which 
1.5 s was steady trim with no information in the data. One cycle 
ot the short-period mode using the batch time-domain parameter 
estimates was 7.5 s. Table 2 contains results at the end of the 14-s 
maneuver from batch time-domain and real-time frequency-domain 
parameter estimations. 

As in the simulated data cases, the standard errors for the pa- 
rameter estimates computed by the real-time parameter estimation 
algorithm were consistent with the accuracy of the parameter esti- 
mates throughout the maneuver. Specifically, the error bounds were 
large at the beginning of the maneuver and smaller as more informa- 
tion was obtained from the data. Except for a few instances that are 
to be expected because of statistical variation, the standard errors 
were representative of the accuracy of the estimated parameters. 

Sampling rate for the flight-test data was 50 Hz, and the recursive 
Fourier transform updates were done at 25 Hz. The same implemen- 
tation and computer were used as before. In this case, the real-time 
parameter estimation algorithm ran roughly 1 4 times faster than real 
time (1 s for a 14-s maneuver) because of the lower rate used for the 
recursive Fourier transform. 

In the final example, control surface inputs measured in flight 
during a piloted longitudinal/lateral tracking task were applied to 
a nonlinear F-16 simulation 11 14 with forward c.g. position (0.2c). 
Figure 0 shows the control surface inputs. The maneuver was ini- 
tiated from a steady trim condition at 10-dcg angle of attack, 
10,000- ft altitude, and Mach 0.32. Simulated output data from the 
nonlinear simulation was corrupted with 20% Gaussian white noise. 
The simulated noisy outputs are plotted in Fig. 10. In this case, 
lateral/directional linear model parameters were estimated in addi- 
tion to the longitudinal model parameters from Eqs. (21) and (22). 


Table 2 FI 8 HARV flight-test results, = 20 deg, 
Au = 24,100 ft, M 0 = 0.34 


Parameter 

Hatch tune domain 
estimate (standard error) 

Rec ill's i v e (re quo ncy -domai n 
estimate (standard error) 

Zu 

-0.218 (0.040) 

-0.200 (0.084) 


1.047 (0.036) 

1.074 (0.052) 


-0.0S7 (0 046) 

-0.041 (0.0S0) 

M a 

-0.649 (0.044) 

-0.509 (0.174) 


-0.063 (0.032) 

-0.177 (0 107) 


-1.257 (0.063) 

-1.415 (0 165) 


F-18 HARV Maneuver 153B 
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Fig. 8 Parameter estimation. 
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Fig. 9 Tracking inputs. 



Time (s) 

Fig. 10 Simulated measured outputs, 20% noise. 



For the lateral/directional aircraft dynamics, the stale vector x and 
input vector w in Eq. ( I ) are defined by 

x = [p p r 4>\ r u = [<$<, (23) 

System matrices containing the model parameters are 

~ Yfi sin a -cos a (g/FJcos# 

_ Lp I'p L r 0 

N fi N n N r 0 

_ 0 l tan$ 0 



Figure 1 1 shows time histones of selected longitudinal and lat- 
eral/directional model parameter estimates, with dashed lines indi- 
cating the ±2a error bounds. Comparison values for the stability 
and control derivatives (plotted as solid lines) were calculated from 
the nonlinear simulation using central linile diflerences with a 1% 
perturbation si/e. Plots lor the model parameters not shown were 
similar. As in the other examples, the performance of the real-time 
estimation algorithm was excellent, in spite of a relatively short 
( 1 0-s) maneuver and low information content longitudinally. Param- 
eter estimates converged to the finite difference values with appro- 
priate standard error estimates. Cycle times for the short-period and 
Dutch roll modes, based on the finite difference parameter values, 
were 3.8 and 2.4 s, respectively. The real-time parameter estimation 
algorithm required 4 s to converge to the finite difference values for 
the model parameters using the tracking inputs measured in flight 
and plotted in Fig. 9. 

The sampling rate for the data was 80 Hz, and the recursive Fourier 
transform updates were done at 40 Hz. For the same implementa- 
tion and computer as before, the real-time estimation algorithm tor 
the longitudinal and lateral/directional models together ran roughly 
eight times faster than real time (1 .2 s for a 1 0-s maneuver). 

Conclusions 

A technique was developed for estimating linear model parame- 
ters in real time using an equation error formulation in the frequency 
domain with a recursive Fourier transform. Simulation and flight- 
test examples demonstrated that the method gives accurate real-time 
estimates of model parameters and standard errors. The examples 
also showed that the advantages of this approach include automatic 



Time (s) 



Fig. 1 1 Parameter estimation. 
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noise filtering, robustness to high noise levels and data dropouts, 
fixed memory requirements regardless of the time record length, no 
bias parameters to estimate, and good performance for low infor- 
mation content iri the data. All of these favorable characten sties, as 
well as low computational requirements, follow from analyzing the 
data in the frequency domain with a recursive Fourier transforma- 
tion using fixed discrete frequencies within the frequency range for 
the dynamic motion of interest. 

The practical applicability of the method was demonstrated us- 
ing a flight-test data example and a nonlinear simulation exam- 
ple using real flight-test tracking inputs implemented by the pi- 
lot. Data information requirements for good parameter estimates 
were found to be low enough that they could be satisfied using 
ordinary pilot inputs measured in flight. The algorithm exhibited 
rapid convergence to accurate parameter values with standard er- 
rors that properly represented the accuracy of the parameter esti- 
mates. Typical convergence times were less than 1 period of the 
dominant dynamic mode. No starting values were required for the 
parameter estimates, no tuning parameters had to be adjusted, and 
there was no temporal or spatial regularization during the parame- 
ter estimation. Parameter estimates and standard errors were based 
on measured data alone. The procedure was shown to have rea- 
sonable computational requirements and ran much faster than real 
time, even when implemented in a high-level language such as 
MATLAB®. 

The technique could be used for dimensional or nondimensional 
parameter estimation and could also be used with general nonlinear 
models, as long as the model is linear in the parameters. Previous 
work 10 has outlined how the technique could be used for real-time 
aerodynamic parameter estimation and flight envelope expansion. 
All states and inputs must be measured, but this should not be a 
problem in modern aircraft with continuous automatic feedback 
control, for which the method is intended. 

The real-time parameter estimation technique studied in this work 
represents a fundamental building block in fulfilling the require- 
ments for adaptive or reconfigurable control. However, each exam- 
ple case included was for a specific flight condition and aircraft 
configuration. Future developments must focus on tracking rapid 
time-varying linear model parameters resulting from changes in the 
flight condition, stores, configuration, remaining fuel, or from var- 
ious types of failures, wear, or damage. In addition, there must be 
some work done to address the question of excitation input design 


when insufficient information content in the data precludes accurate 
parameter estimates. 
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